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Abstract 

We test whether the Tinker & Chen model of Mg II absorption due to the gaseous halo around a galaxy 
' can reproduce absorption in quasar pairs (both lensed and physical) and lensed triples and quads from the 

literature. These quasars exhibit absorption from a total of 38 Mgll systems spanning z = 0.043 — 2.066 
with mean redshift (z) = 1.099 and weighted mean rest-frame equivalent width of 0.87 A. Using the Tinker 
& Chen model to generate simulated sight-lines, we marginalize the unknown parameters of the absorbing 
galaxies: dark matter halo mass, impact parameter, and azimuthal angle on the sky. We determine the 
ability of the model to statistically reproduce the observed variation in Mgll absorption strength between 
paired sight-lines for different values of the gas covering fraction f c and the characteristic length scale £ c , 
within which the variation in absorption equivalent widths between sight-lines exponentially decreases. We 
find a best-fit f c = 0.60 ± 0.15% and l c < 8 /i^q 1 kpc (la confidence limits), with smaller f c allowed at 
larger £ c . At 99.7% confidence, we are able to rule out f c > 0.87 for all values of l c and the region where 
^ ■ l c < 1.0 /ifo 1 kpc and f c < 0.3. 

O ■ 1 Introduction 

Intervening metal absorbers found in the spectra of background quasars have been known to be associated 
with galaxies at intermediate redshifts since the work of Bergeron & Boisse (1991). As a result, quasar sight- 
lines have become a pivotal tool in probing the gaseous halos surrounding these galaxies, using (among other 
transitions) the Mgll AA2796.352, 2803.531 A doublet (e.g., Kacprzak et al. 2011). 

The physical origin of the gaseous halo and its evolving interplay with its host galaxy, however, are not well 
understood. Gas in the halo could be fuel for continued star formation in the galaxy (e.g., Mailer & Bullock 
£ — . 2004), or could be a result of gas blown out by galactic winds (e.g., Steidel et al. 2010), or both. It is clear that 

the gas is clumped and not distributed spherically (e.g., Churchill, Kacprzak & Steidel 2005), but the typical 
^vq , covering fraction is still uncertain. The spatial extent of the Mg n absorbing region has been constrained by 

various studies, although they do not always agree. The extent could be related to the mass of the galaxy 
(e.g., Tinker & Chen 2008, Steidel, Dickinson & Persson 1994, Lovegrove & Simcoe 2011) or to the colour and 
star-formation rate of the galaxy (e.g., Zibetti et al. 2007, Menard et al. 2009 [but see Lopez & Chen 2011], 
Steidel et al. 2010). 

The development of new observational and measurement techniques as well as new models of Mg II absorp- 
tion is required to piece together the picture of the gaseous halo. 

Tinker & Chen (2008) (hereafter TC08) have constructed a semi- analytical model for gas absorbing in the 
Mg II A2796 A transition within a dark matter (DM) halo using an isothermal density profile. The rest-frame 
absorption equivalent width £W(2796) (absorption due to the 2796 A transition of the Mg n doublet at some 
projected distance r from the centre of the absorbing gaseous halo, hereafter EW) is determined by the amount 
of cold gas, in the form of discrete clouds, along the line of sight (LOS) through a halo. They constrained their 
model using a Mg n-LRG cross-correlation function and a Mg II absorbing frequency distribution function. 

In Chen & Tinker (2008), the authors further constrained the TC08 model using a sample of Mgll absorption 
features near known galaxies. The sample consisted of 13 galaxy and absorber pairs and 10 galaxies that do 
not produce Mgll absorption lines to within sensitive upper limits. A maximum likelihood analysis showed 
that the model best describes the observations if the covering fraction is 80-86%. The data had a weighted 
mean rest-frame equivalent width of 0.34 A. 

In this study we present a new test of the TC08 model, using Mg II absorption features in the lines of sight 
to lensed and multiple quasar systems. These systems can provide two or more lines of sight through the same 
galaxy halo. By comparing the coincident and anti-coincident absorption systems in multiple lines of sight, 
constraints on both the covering fraction and coherence scale can be determined (e.g., Ellison et al. 2004). 
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Multiple quasar systems may also be useful in measuring the azimuthal profiles of Mgll absorbers around 
galaxies, a topic studied statistically by Bordoloi et al. (2011). 

We have compiled a non-exhaustive list from the literature of quasar pairs, triples, and quads (hereafter 
referred to as asterisms) where Mgn absorption is observed in one or all lines of sight. Using a statistical 
approach, we ascertain how well the TC08 model can reproduce the observed variation in Mgn absorption 
between sight-lines for various values of the TC08 model parameters. 

We describe our literature sample in § [2] Descriptions of the theoretical framework of the TC08 model and 
the method by which we generate our simulations are given in § EI Our analysis and results are described in § 
HI We discuss our results in § [SJ and summarize our study in § [BJ 

In this paper, we adopt H = 70 km s _1 Mpc -1 , /i7o = H o /70 km s _1 Mpc -1 , Hm = 0.3, = 0.7 (Spergel 
et al. 2007). 

2 Data 

We have tabulated data on the following quasars, for which references are given in Tables Q] and [U 7 
pairs, both physical and lensed, (Q1343+266AB; HE 1104-1805AB; HS 1216+5032AB; Q0957+561AB; SDSS 
J1029+2623BC; SDSS J0904+1512AB; and SDSS J1054+2733AB), 2 triply lensed quasars (Q2237+0305ABC 
and APM08279+5255ABC), and 2 quadruply lensed quasars (H1413+1143ABCD and SDSS J1004+4112ABCD). 
SDSS J1029+2632 B and C are part of a triply lensed source quasar, however, only these two images had avail- 
able spectra and so we treat it as if it were a lensed pair. Q2237+0305 is a quadruply lensed source, but 
only images A, B, and C had available spectra, so we consider it a triple. Also of note in Q2237+0305, the 
absorption system at z = 0.827 only has data for images A and B, therefore we include that one system in the 
pairs data set. 

In the cases of SDSS J1004+4112ABCD, SDSS J1029+2623BC, SDSS J0904+1512AB, and SDSS J1054+2733AB 
we measured the EW ourselves using spectra provided by M. Oguri. For Q2237+0305, we obtained the EW 
by converting the column densities given in Rauch et al. (2002). 

The absorption systems at z = 0.04299 and 0.04431 found in the physical pair HS 1216+5032 have been 
included in this set despite their questionable nature (see § 4.2 in Lopez, Hagen & Reimers 2000). By not 
including the two features in our data set, no major changes to the results were observed. 

Our sample harbours 38 unique absorption features, wherein at least one of the lines of sight in the asterism 
has a rest-frame equivalent width measurement of the Mg n A2796 A transition. The other line(s) of sight in 
each asterism may not have any detectable Mg n absorption, in which case they are assigned a 3<r upper limit 
to their EW (using the RMS noise in the data)Q All data is summarized in Table [T] (pairs) and Table [2] 
(triples and quads) in order of date published (references are found in last column of pairs table, first column 
of triples/quads table). Table Q] lists the object name and components, emission redshift, the lens redshift (if 
applicable), the absorber redshift, the proper separation between the quasar images at the absorber redshift 
(/lyg 1 kpc), and the EW of both lines of sight (measured in A). Table [5] lists similar values, but lists the 
maximum and minimum separation in the object, and does not list the individual proper separations. The 
weighted mean EW in our data set is 0.87 A. 

Our data sample excludes objects with line of sight separations much larger than the scale size of gaseous 
halos described by the TC08 model (proper separations greater than ~300 kpc). 

3 Theoretical Framework and Simulations 

In this analysis we have tested the model developed by Tinker & Chen (2008) using observational data from 
the literature. The TC08 model is a description of how the rest-frame equivalent width of the Mgn A2796 A 
transition changes with respect to the impact parameter between the galaxy and the LOS. Here we describe 
the framework of the TC08 model and how we applied the tabulated data from the literature and coding to 
test it. 

3.1 The TC08 model 

The TC08 model is based on populating a dark matter halo with cold baryonic matter, which is represented 
as Mg II A2796 A absorbing gas. The authors invoke a spherical isothermal density profile to approximate the 

lr The non-detection upper limits were measured in the same manner as the detections in each system. 
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Figure 1: Histogram of all separations in the data set. Bins are 0.25 proper h 70 kpc in the inset and 2.0 proper 
h^Q kpc otherwise. This figure includes all configurations of quads and triples. 



Table 1: Tabulated quasar pairs from the literature that exhibit Mgn absorption at A2796 A. See 
last column for reference. Separations d are quoted in h^ kpc in the proper frame. 



Object 


^abs 


d 


A (A) 


B(A) 


^source i ^lensi R-6f6r6IlC6 


Q1343+266 
9 =975 


0.516 


58.95 


0.65±0.07 


<0.2 


2.03, not a lens, Crotts et al. 1994 


HE 1104-1805 
9 =3'/2 


0.5168 
0.7283 
1.2799 
1.3207 
1.6616 


19.89 

23.24 

9.75 

9.10 

4.76 


0.21±0.02 

0.62±0.02 

<0.05 

0.71±0.02 

0.97±0.05 


<0.23 

<0.27 

0.35±0.06 

0.35±0.10 

0.59±0.11 


2.31, 0.729", Smette ct al. 1995 


HS 1216+5032 
9 =9'/l 


0.04299 a 
0.04431 a 
0.13542 


7.71 
7.93 
21.83 


<0.30 
<0.30 
<0.24 


1.49±0.54 
2.64±0.92 
0.45±0.19 


1.45, not a lens, Lopez et al. 2000 


Q0957+561 
9 =6'/2 


1.3911 


0.235 


2.27±0.01 


2.15±0.01 


1.41, 0.361, Churchill et al. 2003 


SDSS J1029+2632 B-C 
9 =T/768 


0.674 
1.761 
1.910 


10.53 

1.34 

0.81 


1.51±0.07 
0.19±0.04 
0.63±0.10 


1.45±0.24 

<0.4 

<0.6 


2.197, 0.6, Oguri et al. 2008 


J0904+1512 
9 =T/128 


1.2168 
1.6127 
1.6528 


0.52 
0.14 
0.11 


2.0±0.4 
0.3±0.1 
0.3±0.1 


1.3±0.4 

<1.1 

<0.9 


1.826, 0.19, Kayo et al. 2009 


J1054+2733 
9 =T/269 


0.6794 


1.87 


0.77±0.1 


1.07±0.19 


1.452, 0.23, Kayo et al. 2009 



a See the discussion in § 2. 

b The lens redshift was determined by Lidman et al. (2000). 
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Table 2: Tabulated quasar triples and quads from the literature. The value d m ax/d m in indicates the 
largest /smallest proper separations (in units of kpc) between quasar images in the triple or quad at 
the given absorption redshift. For triples the entry is used in the 'D' column. 



Object 


^abs 


A (A) 


B(A) 


C(A) 


D(A) 


dmax 


drain 


H1413+1143 


0.6089 


0.25±0.09 


0.19±0.09 


0.55±0.09 


0.48±0.08 


9.1 


5.1 


Monier et al. 1998 


^source — 2.54 


Zlens — 1-0 












Q2237+0305 


0.566 


0.65±0.07 


0.65±0.06 


0.78±0.07 




0.54 


0.4 




0.827 


0.05±0.06 


0.07±0.01 


b 




0.3 


0.3 


Rauch et al. 2002 


■^■source — 1.69 


z lens = 0.039 












SDSS J1004+4112 


0.676 


0.91±0.03 


<0.2 


<0.2 


1.0±0.1 


101.6 


26.3 




0.726 


<0.1 


<0.2 


<0.3 


2.6±0.1 


93.5 


74.2 




0.749 


<0.1 


<0.1 


<0.2 


1.4±0.1 


89.5 


71.0 




0.833 


<0.1 


0.8±0.1 


<0.2 


<0.4 


77.2 


19.7 




1.022 


<0.1 


0.23±0.03 


<0.1 


<0.2 


52.2 


13.3 




1.083 


<0.1 


<0.2 


<0.2 


1.5±0.1 


44.9 


35.6 




1.226 


<0.1 


<0.1 


<0.2 


0.5±0.1 


31.5 


25.0 




1.258 


<0.1 


<0.1 


<0.2 


0.9±0.1 


28.8 


22.9 


Oguri et al. 2004 


% source — 1.734 


Zlens = 0.68 












APM08279+5255 


1.181 


2.57±0.02 


3.03±0.04 


2.88±0.04 




2.70 


1.07 




1.209 


0.05±0.01 


0.06±0.02 


<0.03 




2.62 


1.03 




1.211 


0.37±0.01 


<0.03 


0.16±0.04 




2.61 


1.03 




1.291 


0.08±0.01 


0.03±0.01 


<0.03 




2.39 


0.95 




1.444 


0.04±0.01 


0.04±0.01 


<0.03 




2.10 


0.80 




1.550 


0.31±0.01 


<0.02 


0.29±0.03 




1.81 


0.72 




1.552 


0.24±0.01 


<0.02 


<0.03 




1.81 


0.71 




1.813 


0.80±0.02 


0.77±0.02 


0.44±0.03 




1.36 


0.54 




2.041 


0.21±0.02 


0.24±0.02 


0.22±0.03 




1.07 


0.42 




2.066 


0.31±0.04 


0.45±0.04 


0.38±0.04 




1.04 


0.41 


Ellison ct al. 2004 


^source = 3.911 


Zlens = 1-062 













a The lensing object is yet to beidentihed, though it is probably at z>l (see Monier et al. 1998). 

b Negligible Mg n absorption was measured in the third image. This feature is used as part of the pair 

sample. 
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mean density of the g function of distance from the centre of the dark matter halo: 



p a (r) = f g G (r 2 + al)- 1 (1) 

where f g is the gas fraction (the fraction of the total baryonic and dark matter mass, Mh, in gas form), ah is 
the core radius, and 

Gq = M h (< R g )/4ir 

R g - a h tan~ 1 (R g /ah) 

is a normalization constant with the arctangent in radians. Only the mass within R g will contribute to the 
density profile. Due to the apparent clumpy nature of the gaseous halo (see, e.g., Rauch, Sargent & Barlow 
1999), the equivalent width of a Mgn A2796 A absorption is directly proportional to the number of absorbing- 
components along the LOS. At radii larger than R g , the density of absorbing clumps drops to zero and no 
absorption is seen. Integrating the individual effects of each clump along a line of sight I through the gaseous 
halo at an impact parameter s relative to the centre of the gas distribution yields 



EW(s\M h ) = EW 



2a c 
M, 



(3) 



EW a a cl f g 2G R l~ s * 

7^+^ tan V 



where a c i and M c \ are the cross section and mean gas mass of an individual clump of gas, EWq is the equivalent 
width per clump, and EW(s\Mh) is the rest frame equivalent width for Mgn A2796 A at s given the dark 
matter halo mass Mh- 

TC08 note that numerical simulations have shown that stable shocks can develop in the centres of dark 
matter halos due to high accretion rates at high halo masses. Above a critical halo mass, the dynamical time- 
scale for matter to accrete onto the dark matter halo is shorter than the post-shock gas cooling time, resulting 
in temperatures at, or above, the virial temperature. The result is a shock-heated region at the centre of the 
gaseous halo in which the cold Mg II absorbing gas can no longer survive. 

To incorporate this effect, TC08 introduce a parameter R s h which is the radius of the transition from 
shock-heated gas to colder gas. At r < R s h the gas has been heated and thus no longer can absorb in the Mgn 
transition effectively. Thus, equation ([3]) is broken into different regions. 

For the cases where R s h = or R s h < s < R g , 



EW(s\M h ) = ^^t,n-\ (4) 

If the impact parameter is within the shock- heated region (s < R s h), then the effect of the reduced Mgn 
absorbing gas is subtracted from the total: 



2A EW G x \ R 2 sh - s 2 



EW(s\M h ) = EW [Rsh=0) - f hot J^ tan- 1 J f (5) 



where fhot — 1 ~ fcoid and f co id is the fraction of gas within the shock-heated region that is still able to absorb 
in Mgn (cold gas). If the entire gaseous halo has been shock- heated (R s h — R g ) then we measure the EW by 



2G 



EW(s\M h ) = f cold A EW ^=JL= tan' 1 J J ; \ (6) 



as only a fraction f co id of the gas within the halo will absorb at the Mgn doublet transitions. Of course, if 
s > R g then EW(s\M h ) = 0. 

In the above formulae TC08 grouped the parameters EWq, <r c i,f g , and M d into the single parameter, A EW , 
which represents the mean absorption equivalent width per unit total surface mass density of the cold gas: 

EW 

Aew = tz i ( • (7) 



■5 




Figure 2: The equivalent width of the Mgn 2796A transition as a function of impact parameter (in comoving 
reference frame) predicted by the TC08 model. Each curve represents a different dark matter halo mass. The 
legend provides the curve for each of the values for mass, measured in solar masses and in the log. In the top 
figure, the effects of the shock heated region are manifested as a reduction in EW at low impact parameters 
for masses between log M h /M e = 11.00 and log M h /M Q = 11.75. At log M h /M Q = 12.00, the entire gaseous 
halo becomes shock heated and the underlying form of the density profile is recovered, although, at reduced 
EW. In the bottom figure, the model curve for the mass log M h /M Q — 11.75 is plotted again for reference. 
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Figure [2] shows the EW(s\Mh) profile for various values of mass. Note the cavity created at small impact 
parameters for halos with shock-heated regions. 

The value of R g (as well as ah and R s h) are dependent on the mass of the halo. Following TC08 we set 

R g = 8O(M /l /lO 12 M ) 1/3 /i~ 1 kpc (comoving) (8) 

where h = Hq/(100 km s _1 Mpc -1 ) and ah = 0.2R g (TC08 note that the exact choice of ah has little effect on 
their results). Converting to the proper frame and using our Hubble constant, the radius of the gaseous halo is 

where = Hq/(70 km s^ 1 Mpc -1 ) and z a bs is the redshift of the halo. Hereafter, all coordinates will be in 
the proper frame. 

TC08 have parameterized the value R s h as 

^ = R° sh + Ish Iog(il4/1O 12 M ) (10) 

Rg 

and forced the bounds < R s h/R g < 1. Any values of R s h/R g < correspond to a gaseous halo that has 
no shock- heated region (i.e., R s h = 0), and R s h/R g > 1 indicates a gaseous halo that is entirely shock- heated 
(i.e., R sh = R g ). 

TC08 utilized a Monte Carlo Markov Chain technique to determine the best fit values to the 4 free pa- 
rameters in the model. We adopt the same parameters. These are Aew = cm" g - 1 = 374 hro A 
cm 2 g _1 , f co id — 0.061, R® h = 1.02 and 7^ = 1.03. With these parameters, a shock radius of R s h = has 
mass logM^/M® < 11.0 and R s h = Rg has mass lagMh/M® > 12.0. The mass range over which the shock 
heated region appears and encompasses the entire gaseous halo is small in comparison to the mass range used 
in TC08's analysis. 

The TC08 model was tested against 13 galaxy/absorber pairs in Chen & Tinker (2008) where a cosmic 
scatter of 0.25 in the log was observed. With such a large cosmic scatter, a pair of sight-lines with increasingly 
smaller separation still has the potential to have significantly different absorption equivalent widths. We expect 
there to be some scale at small separations over which the absorption equivalent width will not vary between 
sight-lines. To account for this we introduce a characteristic length, denoted £ c . Below £ c , the cosmic scatter 
is suppressed exponentially so that at zero sight-line separation, the difference in cosmic scatter between sight- 
lines is zero (see section 13.2.21 for more details) . 

We include a halo absorbing fraction fh{Mh) as defined in TC08. This is the probability that a dark matter 
halo has any gaseous halo capable of absorbing at all. More details on this parameter can be found in section 

We also introduce a covering fraction term, f c , which measures the fraction of sight-lines within R g that 
do not show absorption (see section 13". 2. 41 for more details). 

There are a total of seven free parameters required to fully describe a gaseous halo in our model. TC08 
determined the values for five of these parameters. In this study we have varied the remaining two, the 
characteristic length and covering fraction, in order to find the best fit values. The characteristic length £ c is 
varied over 0.1-100 /iff, 1 kpc and the covering fraction f c is varied over 0.01 — 1.0. 

3.2 Simulating Paired Quasar Sight-Lines 

We aim to determine if the variation in EW between sight-lines in the observed data set (Tables [1] and ^ 
is consistent with the TC08 model, and over what range of characteristic length and covering fraction that 
consistency holds by simulating sight-lines in a gaseous halo described by the TC08 model. 

We begin by choosing a mass at random from a halo mass function. The mass defines the radius of the 
gaseous halo, and the size of the shock heated region. We take a random impact parameter from within the 
gaseous halo and calculate the equivalent width predicted by the TC08 model, based on equations (gJIHl). We 
search through all EW values measured to better than 3a in the observed data set for a match to the predicted 
EW0 The matched EW represents one LOS from the corresponding quasar asterism. We use the TC08 model 
to then predict the EWs of the other lines of sight in the asterism, given the mass of the halo and the impact 

2 A single predicted EW from the TC08 model can be matched to multiple observed EWs. 



7 



parameter that was used in the original match. For each line of sight in the observed data set with an EW 
measured to better than 3c, we compile 900 EWs that match to within observational error. For each of the 900 
matches, we then use the TC08 model to predict the EWs for each of the other lines of sight in the asterism, 
given an assumed characteristic length scale t c and covering fraction f c of the absorption. Using the resulting 
distribution of predicted EWs, we calculate the probability that the TC08 model can accurately predict the 
observed variation between lines of sight for each combination of £ c and f c . 

This statistical approach marginalizes the mass of the dark matter halo, the impact parameter of the sight- 
line relative to the centre of the gaseous halo, and the orientation of the asterism within the halo. Below we 
describe how we determine these parameters, and how we generate simulated samples of EWs to compare with 
the observed EWs. 



3.2.1 Selecting Mass 

We select a mass randomly from a weighted halo mass function. The weight we have already defined as fh(Mh), 
the probability that a given dark matter halo of mass Mh will have a gaseous halo from which absorption can 
be measured, as not all halos contain a gaseous region. The halo mass function quantifies the number of 
dark matter halos as a function of halo mass. TC08 utilize the halo mass function of Warren et al. (2006), 
determined by numerical simulations of structure growth and halo formation. Multiplying the mass function 
by fh results in an absorption weighted halo mass function, allowing proper sampling of the mass range for 
dark matter halos that would be observed in our data set. 

In TC08, however, //, is degenerate with the covering fraction of the gaseous halo. So the authors combine 
the probability that a given dark matter halo will have a gaseous halo that will absorb in Mgn, and the 
probability that a given line of sight inside the gaseous region will produce Mgn absorption. In this paper, 
we are required to explain cases of absorption in only one of two lines of sight that lie within R g . Since we 
require there to already be a match to one line of sight in an asterism (discussed in § I3.2.2|) . we are not double 
counting the effects of covering fraction by using fh and f c . 

Instead of parameterizing the mass dependence of fh, TC08 specified the values of fh at four different 
masses and interpolated between them. The four masses used were log Mi/M Q = 10.0, 11.33, 12.66, and 
14.0. The four corresponding values of fh are denoted /i,/2,/3 and f±. In TC08, the best fit values for fh 
at each log Mi were determined; their values were: log/i = —9.27, log/2 = —0.205, log/3 = —0.006, and 



log fi = —0.168. We adopt the same parameters and interpolate linearly between them; the plot in Figure 3(a 
shows the dependence of fh on mass (both in the log). 

For the sake of simplicity, we approximate the halo mass function as a power law (dn/dMh — Mr a ) that we 
sample over the range 10.75 < log M^/M < 14.75. In Warren et al. (2006), the parameterization is actually 
a power law combined with an exponential decay at high masses (> 1O 145 /i _1 M0), however, our sampling 
of masses in this range is negligible and so we consider our approximation to be sufficiently accurate. From 
the power law portion of the Warren et al. (2006) halo mass function, we measure a — 1.88. Therefore, our 
weighted halo mass function is defined as 

dn w = fhMh a (u) 
dM Sl^thM^dMh 

such that we have a normalized weighted halo mass function with sampling range 10.75 < logM/j/M© < 14.75. 
A plot of dn w /dMh is located in Figure [3 (b)| where n w represents the weighted number of absorbers. 

As noted in TC08, the weighted mass function peaks at log M^/M© ~ 12, and nearly all halos of 11.5 < 
log Mh/M® < 12.5 are expected to host Mgn absorption. At lower masses, the halo mass function has a 
relatively high normalization but the low value of fh makes the probability of finding an Mg 11 absorber in the 
dark matter halo low, and thus we have a small chance of observing them. At higher masses, fh is relatively 
high but the number of halos is much lower, making the probability of observing a Mg 11 absorber at such halo 
masses very unlikely. 

The Cumulative Distribution Function (CDF) of the weighted halo mass function is plotted in Figure [4] 
This plot shows that the masses we are most likely to be studying are approximately 11 < log Mh/M Q < 13. 
All other masses have a low probability of exhibiting Mgn absorption, although we still sample over the range 
10.75 < logM/j/MQ4 < 14.75. The CDF allows us to randomly select dark matter halo masses to properly 
represent the statistical distribution of dark matter halos in the universe that will harbour Mgn absorbing 
gas. This model assumes no evolution of the halo mass distribution of absorbers with redshift, which has been 
suggested by observational evidence in Lundgren et al. (2011). 
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Figure 3: Top The probability that a dark matter halo of mass Mh will show Mgn absorption, fh, plotted 
against halo mass. Note that at low masses, the probability is very low. Much of the Mgll absorbing gas is 
found in dark matter halos with masses > 10 11 M Q . The form of fh was found by linearly interpolating between 
4 values determined by the best fitting procedure in TC08. Bottom The weighted halo mass function plotted 
against halo mass. The low value of fh at low masses forces dn w /dMh to be relatively low at low masses, and 
the power law form of the halo mass function forces dn w /dMh to be relatively low at high masses. 
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Figure 4: The Cumulative Distribution Function found by integrating over dn w /dMh for all values of mass. 
Masses in the range 11 < log M^/M© < 13 are much more likely to exhibit Mgn absorption compared to all 
other masses (high or low). 
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3.2.2 Simulations for Pairs 



The location (i.e., impact parameter) and orientation of an observed quasar asterism in reference to its absorbing 
gaseous halo is unknown. Here we describe how we generate a random distribution of simulated data within the 
constraints of the TC08 model that incorporates our ignorance of the location of the Mgn-absorbing galaxy. 
We also describe how we test this distribution of simulated data against the observed data in Table [TJ It is 
important to note that when simulating EWs we search for matches to both LOSA and LOSB in the observed 
quasar pair. The only stipulation is that the observed EW in the LOS must be measured to better than 3a 
(i.e., we do not match to upper-limits). 

Each run of simulations begins by choosing and fixing the characteristic length and covering fraction com- 
bination to be used during the run. These are chosen from within the ranges given at the end of § 13.11 We 
then generate a halo mass Mh by selecting at random from the weighted halo mass function described in the 
previous section. The mass is used to calculate R g using equation (|9]), R s h using equation (flO| . and Go using 
equation (|2|). An impact parameter s a is then randomly generated within the cloud of radius R g , defining a 
location a (all parameters defined at this location will have the subscript a). Using these model parameters, 
impact parameter s a , and equations (gHH]) we generate EW(s a \Mh), the EW due to Mgn absorption at impact 
parameter s a given the halo mass Mh- 

In Chen & Tinker (2008), the authors found the model used in this study had a cosmic scatter (a cs ) of 0.25 
dex. To incorporate a cs , we must scatter the simulated EW(s a \Mf l ) individually for each absorption system 
in our quasar sample. For each line of sight in a pair that is measured to better than 3cr the scattered EW is 

log£W a = \ogEW(s a \M h ) + a a (l - e"^) + a x e~ d '^ (12) 

where a a is the cosmic scatter at the location a, x denotes the point on the sky midway between the lines of 
sight of the observed pair, a x is the cosmic scatter at x, and d x is the proper separation between a and x at the 
absorber redshift. The values a a and a x are chosen at random from a Gaussian with ji = and a — 0.25 dex. 
The characteristic length l c is one of the variable parameters in the model. Note that when d x ^> i c , the a x 
term is negligible and the cosmic scatter is a a at a. If d x £ c then a x dominates the cosmic scatter at a, while 
the a a term is negligible. This relationship is used because we assume there is some scale on which cosmic 
scatter should not apply, and the two lines of sight will be, effectively, identical. However, if the separation 
of the pair is much larger than the characteristic length, then cosmic scatter may be applied to each LOS 
calculation individually. 

Note that the value d x in equation (|12p . is unique to each absorption system, and thus the value of EW a 
is calculated individually for each entry in Table [T] before checking to see if it is a match. 

Using EW(s a \Mh) and equation (p~2|) . we search the different lines of sight in the pair sample for an EW that 
matches the scattered EW a (note that we can match either LOSA or LOSB to this value). To be considered 
a match the simulated value must be within ±\/3<7 of the given observed EW (where a is the observational 
uncertainty on the EW; see Table H} |f| 

In cases where the EW in both lines of sight at a given redshift in an observed pair are measured at > 3<r, 
we check for matches with both the first line of sight (LOSA) and the second line of sight (LOSB) in the 
pair individually (see Table [1] for which LOS is 'A'). In cases where only the EW in LOSA is measured to a 
statistically significant amount, and LOSB has been quoted as an upper limit, we search for matches to LOSA 
only (and vice versa). The cases for triplets and quads will be described later (see section l3.2.3[) . It is entirely 
possible that a randomly generated EW(s a \Mh) does not scatter to match the EW at either LOSA or LOSB 
in a given absorption system. 

When a match to an observed LOS is made, we use the TC08 model to compute a second EW, located at 
s b , which is a distance d away from the matched LOS located at s a (The distance d is the separation listed in 
column 3 of Table [T]) . This second EW is generated by 

\ogEW b = \ogEW(s b \M h ) + a b {l - e'^) + a x e~ d ^^. (13) 

The impact parameter s b is the location of the second line of sight in the quasar pair, and is found by placing 
LOSB at a separation of d from LOSA and randomly rotating it in the plane of the sky an angle of < <fi < 2ir. 
The value EW(s b \Mh) is then determined using equations ((310 • The value a b is found in the same way as a a 
and The resulting EW b value is the TC08 predicted EW for a LOSB in a quasar pair, given a match to 
LOSA (EW a ) and the proper separation of the pair. Again, note that we also allow for the reverse: matching 

3 This range is used because a box function of width 2<t\/3 has the same variance as a Gaussian function with variance a 2 . 
4 The value of a x is the same for both EW a and EWi,. 
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to LOSB and predicting LOSA. In the limit of d x = 0, s a = Sb (a and b represent the same point) so that 
EW(s a \M h ) = EW{s b \M h ) and \ogEW a = \ogEW(s a \M h ) + a x = log£W&, as expected. Application of the 
covering fraction is discussed in § 13.2.41 

For each line of sight in the 18 absorption systems in the quasar pair sample that was not measured as 
an upper limit, we tabulated 900 simulated (EW a ,EWb) pairs. Each set of 900 contains simulated EWs that 
have the same corresponding proper separation as the observed absorption system, but random halo masses, 
projected impact parameters, and orientations on the plane of the sky. Once each LOS in the pairs sample 
has tabulated 900 pairs, we re-run the simulations for a new combination of characteristic length and covering 
fraction. 

3.2.3 Simulations for Triples and Quads 

Applying the model to a quasar triple or quad requires a slightly different approach, in that there are now 
three or four lines of sight to consider when predicting EWs. 

The equation for EW a is conceptually the same as for pairs. A random location in the absorbing halo is 
chosen. We define its impact parameter to be s a , and we have: 

log EW a = \ogEW{s a \M h ) + <7 Q (1 - e -<W'«) + a x e- d ^ e ° (14) 

where a a is the cosmic scatter at the location given by s a , d ax is the proper distance from the LOS to the 
projected centre x of the quasar asterism (see Tabled, and <r x is the cosmic scatter at x. The values a a and 
o~ x are chosen at random from a Gaussian with /i = and a — 0.25 dex. 

In each triple/quad there are multiple lines of sight per absorption system. We check each LOS for a match 
to the random EW, with the same stipulations for a match as used by the pairs simulation. When a match 
is made to one of the observed lines of sight in the asterism, we place a scaled, rotated copy of that asterism 
in the gas cloud, with the LOS to which a match has been made located at impact parameter s a . The scaling 
takes into account that the proper distances between the lines of sight are dependent on the redshift of the 
absorption feature. The asterism copy is rotated by a randomly chosen angle <f> between and 2-7T. To calculate 
the EW n values of the other lines of sight in the asterism, we use 

log EW n = log EW{s n \M) + o»(l - e -<W4 ) + , ie -4./«. (15) 

where n = b, c, or d, d nx is the equivalent of d ax for each sight-line and o~ n is the equivalent of a a for each 
sight-line. The values EW(s n \M n ) are determined by equations (J3HHI) - The fourth EW {EWa) is only generated 
if the object is a quad. This method generates three or four EWs in the same geometry as an observed object, 
at a random point inside a gaseous halo described by one randomly selected mass. For each LOS with an EW 
measured at > 3a in each object, we tabulate 900 sets of simulated EWs. 

As a sanity check, in Figure [5] we have plotted the distribution of masses generated by the simulations 
that ended in a match to the EW in a line of sight from the data set. We plotted the distribution for two 
different lines of sight, a weak absorber and a strong absorber^ Both absorbers were from the APM08279+5255 
simulations. 

The two mass distributions are clearly distinct, but keeping Figure [2] in mind, they make sense. The less 
massive halos have TC08 modeled EW profiles that are more able to achieve larger EW values. In randomly 
generating masses and looking for matches, it makes sense that the distribution for a strong absorber will 
include more massive halos. 

3.2.4 Covering Fraction 

We have not yet accounted for the covering fraction. Effectively, the above description is one of simulating line 
of sight EWs with 100% covering fraction. As discussed in section Q] the measurements of f c have a wide range 
of values. We test a range of covering fractions from 0.01 to 1.0, which encompass the results of Kacprzak 
et al. (2008) (/ c = 0.5 for 0.02 A < EW < 3.0 A) and Barton & Cooke (2009) (f c < 0.4 for weighted mean 
rest-frame EW — 1.29 A) as well as the covering fraction measured from this model by Chen & Tinker (2008), 
f c = 0.85 for weighted mean rest-frame EW = 0.34 A. Note that the covering fraction measured in these 
studies depends on the EW range of the data set, with smaller f c at higher EW. 

Each simulation run has a fixed characteristic length and covering fraction. When a match to an observed 
sight-line is made by scattering EW(s a \Mh), each of the other sight-lines in the quasar asterism would then 

5 The dividing line between 'weak' and 'strong' has been defined in the literature as EW = 0.3A. 
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Figure 5: The distribution of masses generated by our simulations that resulted in a match to a line of sight. 
The white histogram is for a strong absorber with EW = 0.80 ± 0.2 A. The black histogram is for a weak 
absorber with EW = 0.05 ± 0.01 A. Both absorbers are from APM08279+5255 (Ellison et al. 2004). Each 
histogram totals 900 different masses. 
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have a f c chance of actually absorbing in Mg II. (We do not apply the covering fraction to sight-line A because 
we are measuring the ability of the TC08 model to predict the EW in other sight-lines given the matched EW 
at LOSA.) 

For pairs, to account for the covering fraction in the LOS located at position b, we randomly generate a 
number between and 1 for each of the 900 (EW a , EWb) simulated pairs. If the random number is above the 
covering fraction value, we reset EWb to the value EWl given by 

EW' h = [EW a + a a ]e~ d ^. (16) 

where EW a is the matched EW, a a is the cosmic scatter at a, and d = 2d x is the separation between LOSA 
and LOSB. In the case of triples and quads, we also recalculate the EW at each of the lines of sight predicted 
by the TC08 model using a random number approach, however, we do this relative to the EW at the centre of 
the asterism: 

EW' n = [EW(s x \M h ) + o x \ e -<W*c (17) 

where n = b, c or d, EW(s x \Mf l ) is the EW at x, a x is the cosmic scatter at x, and d nx is the distance from 
LOSn to x. 

As a result, for all quasar asterisms in our tabulated data from the literature we have created 900 unique 
and independent (EW a , EW n ) sets at each value of £ c and f c with the value EW a matched to one of the lines 
of sight of the parent quasar asterism. Each of the 900 simulated sets was generated for a randomly chosen 
dark matter halo mass, projected impact parameter, orientation, and probability f c of absorbing in Mgnat a 
given location. 



4 Analysis: Comparing the Simulations with the Observations 

In this section we describe how we analyze the simulations. We use the Kolmogorov-Smirnov (KS) test (in § 4.3) 
to compare the simulations to what is expected if the TC08 model exactly predicts the observed distribution of 
equivalent widths. To that end, we take into account the differences in relative size of the gaseous halos used 
to generate the simulated EWs. We use the cumulative distribution function of the fractional probabilities for 
each paired sight-line (§ 4.2) as a measurement of the ability of the TC08 model to predict the observed EWs. 

The original TC08 model was constrained using a sample ranging in redshift over 0.2067 < z < 0.892. We 
provide the KS test for both our entire sample, and a redshift limited sample in the same range originally used 
in the TC08 model. 

For illustrative purposes, below we denote an observed EW as EW ±a and its corresponding 900 simulated 
EWs as EWf m . 

4.1 Relative Size Considerations 

Because gaseous halos of different masses and thus sizes were used to create the distribution of EWp m for each 
EW Q ± cr , we must take into account the relative probability each individual halo has of producing absorption 
that matches an observation (call it P'j). Each halo mass has a unique gaseous halo size as determined by 
R g , a unique equivalent width profile (see Figure [2]), and a cosmic scatter of 0.25 dex in the EW at any given 
position. Therefore, each impact parameter in each gaseous halo has a different probability of producing an 
absorption feature that matches an absorption feature from the observed data set. 

For any one impact parameter s relative to the centre of the jth gaseous halo of mass Mj and radius Rj , 
we weight that impact parameter following 

2ns ■P abs {EW ±o \s) (18) 

where P a bs{EW ± a \s) is the probability of generating an equivalent width match to an observation of 
EW ± cr , given s and the equivalent width profile of the gaseous halo. 

To determine P a bs{EW ±a \s), we first need P a bs(s), the probability distribution of EWs at s for a given 
Mj in the TC08 model. The exact form of P a bs(s) is 

p M n f -[logEW - logEW r {s\M 3 )f \ 

Pabs{s) = £>exp I — I (19) 

where EW r (s\Mj) is the simulated equivalent width at the given impact parameter (see equations 0HH]), a cs — 
0.25 and the constant D is a normalization factor. We approximate the Gaussian in equation (j!9[) as a 
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Figure 6: The dotted line labeled 'EW' represents the unscattered TC08 model for M = 10 n - 2 Af Q . The solid 
line labeled 'EW max .' represents the ler cosmic scatter above the TC08 model; the cosmic scatter below the 
TC08 model is also whosn as a solid line. An observed Mgn absorber is labeled EW ±1.73a, where -\/3 ~ 1.73. 
The region where the the observed Mg n system and the scattered TC08 model overlap represents the relative 
probability that a gaseous halo defined by the above mass will generate a simulated equivalent width that 
matches the observed value. 



uniform probability distribution in log r between log EW m i n — log EW r {s\M) — cr C sV3 and log EWmax = 
log EW r (s\M) + a cs \/3. To determine the relative probability that the individual impact parameter s will 
produce an absorption feature that matches an observation EW Q ± a we measure the area between EW m i n 
and EWmax that overlaps the allowed range around the observation. For each individual impact parameter, 
the area where the bounds of the observed EW D ± a a and the bounds of the TC08 model overlap represents 
the relative probability of producing an absorption feature that would be considered a match: 



P abs (EW ±a \s) 



dlogEW 



(20) 



log EW m ax ~ log EWmin 

where a = max(log EW m i n , log EW — cr v3) 

and b = min(logP W max , log EW a + a a V3). 

In Figure [5] we show an example of this calculation. An observed equivalent width, labeled EW Q , has some 
range over which we consider it to be matched, determined by the observed uncertainty ±cr\/3- In the case of 
Figure [6j the equivalent width profile EW was created using a halo mass of log Mh/M® = 11.5. 

As Figure |6] shows, it is possible to match EW over a range of impact parameters, each of which contributes 
to the probability of matching the observation in a halo of size Rf, we thus sum over all impact parameters 
within the jth gaseous halo to find the relative probability of obtaining a match to the observed value EW 
in that halo. We then divide by the sum of the relative probabilities of all N simulated matches to find P', 
which is the relative probability of matching the observed EW D in a halo with radius Rj : 



P' 

3 



N / 

}fracY / ZksSs ■ P abs (EW ± a \s)J2 

s<Rj fe=l \s<R k 



2tts Ss-P abs (EW ±a \s) 



(21) 
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where N is the number of simulated matches. 

Each of the 900 EWp m are no longer weighted equally at 1/900, but now weighted as Pj/900. As a simple 
example, if P a i, s (EW ± cr |s) = fixed independent of s or Mj then the weight of each halo is just cx R^, and 
the relative weight of each halo is Pj = P|/Q^j R^)- 

4.2 Fractional Probability Curves 

For each choice of £ c and f c , we measure how well the simulated EWs match the observed measurements by 
looking at the cumulative distribution function (CDF) of the fractional probabilities for each paired sight-line. 
Conceptually, the fractional probability is the fraction of the 900 simulated EWs which are equal to or less 
the actual measured EW of LOSB in that paired sight-line. If the model being tested matches the 
observations well, then the fractional probabilities generated by that model for our sample of paired sight-lines 
should be uniformly distributed between and 1. 

However, some of the EW values are measured as upper limits and the remaining have uncertainty values. 
Measuring simply the number of values below EW does not include the observational uncertainties associated 
with it. To account for this, we compared the weighted cumulative distribution of the EWp m (weighted by 
Pj)to the cumulative distribution of a Gaussian with fi = EW and a = <x,-; i.e., an error function erf(/i, er)0 
The erf (/i = EW , a — a ) values plotted as a function of the weighted cumulative distribution of the EWp m 
provide a fractional probability curve for each paired sight-line. That curve combines the fractional probability 
at each equivalent width with the relative probability of each equivalent width being the true value given the 
observed EW a and <j . Thus, instead of generating the CDF of fractional probabilities for all sight-lines by 
summing step functions corresponding to a single fractional probability for each paired sight-line, we generate 
it by summing the fractional probability curves for each paired sight-line. 

4.3 The Kolmogorov-Smirnov Test 

We can use the Kolmogorov-Smirnov (Press et al. 2007) test to see how well the simulated TC08 model 
equivalent widths match the observations. 

The Kolmogorov-Smirnov (KS) test is a way of measuring the difference between two CDFs. The amount 
by which the CDF Sn(x) of some given data set differs from a known distribution denoted P(x) is measured 
by finding the maximum value of the absolute difference between the two CDFs. This is denoted D and defined 
as 

D = max \S N (x) - P(x)\ . (22) 

To avoid decreased sensitivity near P(x) — or 1, we use the variant of the KS test known as the Kuiper 
statistic (Press et al. 2007). The Kuiper statistic V measures the maximum value of the difference in the CDFs 
above and below P(x): 

V = max [S N (x) - P(x)} + max [P(x) - S N (x)} . (23) 

If the simulated data generated by the TC08 model above matches well the observed data points from the 
literature, we would expect the sum of the fractional probability curves of the individual observed EWs to 
create a cumulative distribution function that is similar to P(x) — x. We can measure the deviation from 
P(x) — x using the Kuiper statistic. The model parameters that yield the smallest value of V therefore 
represent the parameters that best match the observed values of EW. 

We seek to determine how well each combination of the model parameters £ c and f c reproduces the observed 
data. For this analysis, we separated the triples and quads into sets of pairs. A triple quasar can be considered 
a set of 6 individual quasar pairs. A quad can be considered a set of 12 individual quasar pairsH In doing so 
we do not lose information regarding the two dimensionality of the asterisms, as the EWs of each line of sight 
were generated keeping the form of the asterism in the gaseous cloud consistent. We only break the asterisms 
apart in order to calculate individual fractional probability curves. 

In total, there are 25 fractional probability curves from the pair sample, 48 from the triples sample, and 39 
fractional probability curves from the quads sample, with a total set of 112 individual fractional probability 
curves per combination of (£ c , f c ). 

6 The fraction equal to or greater than the actual measured EW is an equally good statistic. 

7 If the observed EW is actually an upper limit, we set fi = and a = <r , where a is equal to one third the value of the upper 
limit quoted in Table [l] 

8 This approach is conservative in the sense that we do not assume that matching LOSA and predicting LOSB yields the same 
results as matching LOSB and predicting LOSA. 



10 




0.0 Lil_, , , I , , , I , , , I , , , I , , , L 

0.0 0.2 0.4 0.6 0.8 1.0 

fractional probability 

Figure 7: An example of a bad fit to the null hypothesis. The fractional probability curve is for t c = 10.0/ifokpc 
and f c = 0.9. The p- value probability is effectively 0.0 indicating that, with these model parameters, the TC08 
model cannot reproduce the observations. 
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Figure 8: The fractional probability curve for l c = 2.0 /ifo 1 kpc and f c = 0.2. This is a much better fit than the 
model parameters in Figure [71 and the p-value probability of ~ 0.13 indicates that we cannot rule out these 
parameters to better than 3<r. 
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Figure 10: Contour plot showing the results of applying the KS test with Kuiper variant to each combination 
of (£ c , f c ). The la, 2a, and 3tr contours represent regions outside of which we can rule out the TC08 model to 
68.3%, 95.4%, and 99.7% confidence, respectively. The black square denotes the best-fit characteristic length 
£ c = 0.5 /ifg 1 kpc and covering fraction f c = 0.6. 

For Figures [JJ [51 and [HI the dotted line represents the null hypothesis: if the model exactly matches the 
observed data set, we expect the fractional probabilities to create a curve of P(x) = x. The solid line is the 
fractional probability curve from the simulated data. The vertical dashed lines show the positions of the largest 
positive and negative deviations of the histogram from the P{x) = x curve. The V statistic (see equation 123 j) 
and the p- value are both labeled in the top left corner of the graphs. The p- value is the probability, given the 
data, of determining a value for V greater than the one found, adopting the null hypothesis. 

For each point in (£ c , f c ) space we calculate the p-value. A contour plot in (i c , f c ) space is shown in Fig. 
[TOl We are able to reject the TC08 model in all parameter space outside of the 3a contour. The 2a and la 
lines are also plotted. The best fit to the observed data occurs at £ c = 0.5 kpc and f c = 0.6, measured for 
a mean rest-frame EW = 0.87A. 

4.4 A Redshift Limited Sample 

In TC08, the model was initially constrained using both a Mgn-LRG cross-correlation function built from a 
sample that ranged from z — 0.35 — 0.80, and a Mgn absorption frequency distribution function compiled from 
two sources that spanned the range z = 0.2 — 2.2 (see § 3 of Tinker & Chen 2008). The TC08 model was then 
tested using a sample of 13 galaxy /absorber pairs with an average absorber redshift of (z a b s ) = 0.3818 ranging 
z = 0.2067 — 0.892 (Chen & Tinker 2008). The study yielded a cosmic scatter of 0.25 dex, which we have 
adopted in this study. 

The mean redshift of the Mg n absorbers in our sample was (z a bs) = 1.099 with a range of z = 0.043 — 2.066. 
Our sample probes to smaller redshifts than the TC08 model was initially constrained with, but the contribution 
of our sample to the z < 0.2 regime is only 2 out of 38 absorbers or ~5% of the sample. We were interested 
to see if using the subset of Mg II absorbers in our sample that were within the redshift range used to observe 
the 0.25 dex cosmic scatter would change the constraints on the model parameters we varied. 
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Figure 11: Re-plotting the contours using only Mgil absorbers from our sample that lie within the redshift 
range 0.2067 < z < 0.892. This range matches that used by Chen & Tinker (2008) when testing the TC08 
model. The black square denotes the best-fit characteristic length t c — 1.5 /lyg 1 kpc and covering fraction 
fc = 0.6. 



In Figure QTJ we re-calculated the p- value probabilities and then re-plotted the contours using only simula- 
tions generated from observations whose Mgn absorption feature lies within 0.2067 < z < 0.892. By limiting 
the redshift range, we reduced our absorber sample from 38 to 12. 

The redshift limited sample is much less constraining than the entire sample (Figure ITUjl . This is most 
likely due to the drastically reduced sample size. The KS test with Kuipcr variant is heavily dependent on the 
number of data points used, and the lack of constraint in Figure QTJ is attributed to the lack of data points 
involved. 



5 Discussion 

Our analysis has ruled out significant portions of the parameter space we tested. All characteristic lengths 
larger than 20 h^Q kpc are ruled out, regardless of covering fraction. This result indicates that the underlying 
form of the TC08 model does not represent the distribution of the Mgn absorbing gas. When the characteristic 
length is so large, it forces both sight-lines in a pair to be scattered by the same value a x . This preserves the 
underlying form of the TC08 model (i.e., the 'unscattered' form), and is thus 'halo-to-halo' scatter; all sight- 
lines in one gaseous halo are scattered the same amount. This is different from 'point-to-point' scatter, where 
all sight-lines in one gaseous halo are scattered by different amounts. By ruling out large characteristic lengths, 
we have found that the TC08 model requires the cosmic scatter to be applied point-to-point, and thus the 
underlying analytical form of the model cannot fully reproduce the observed variation in Mgn absorption 
between sight-lines. 

Given the above argument, it is not surprising that we have also ruled out covering fractions above 0.87 
(for all characteristic lengths). At covering fractions above 0.87, the probability that any sight-line will be set 
to zero is low, and thus we are forcing the use of the underlying analytical form of the TC08 model (see Figure 
[2]) . Observations show that not all lines of sight below the estimated R g show Mg II absorption, and thus it is 
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not surprising that high covering fractions are ruled out. 

There is also a region at low characteristic lengths and low covering fractions that has been ruled out (see 
the bottom left of Figure IT0| . At covering fractions f c < 0.3, only the matched sight-line is guaranteed to have 
absorption; most other sight-lines only have residual absorption governed by the characteristic length, as set 
out by equations [TH] and [T7J If the separation between sight-lines is much larger than the characteristic length 
of the simulation, most b, c, or d sight-lines will be set to zero. In the literature this is clearly not the case, 
and thus we are able to rule out low covering fractions when the characteristic length is much smaller than the 
average separation of the data set. 

As the characteristic length increases, the exponential factors of equations 1161 and 1 1 71 are no longer approx- 
imately zero. As a result, even though most b, c, or d sight-lines are randomly chosen to be set to zero (at 
f c < 0.3), the large characteristic length forces many to be non-zero. This effect is clearly visible in Figure HU1 
where we can no longer rule out low covering fractions between 2.0 < £ c /hm kpc < 20. 

6 Summary 

We have tested whether the model due to Tinker & Chen of Mg II absorption from the gaseous halo around a 
galaxy can reproduce absorption in quasar pairs (both lensed and physical) and lensed triples and quads from 
the literature. Specifically, we have determined what the acceptable values of the characteristic length and 
covering fraction. 

The TC08 model predicts the magnitude of Mg n absorption due to a gaseous halo due to a galaxy as a 
function of the impact parameter of the incident sight-line with respect to the centre of the gaseous halo. In 
this paper, we have tested the TC08 model using a subset of quasar pairs (both lensed and physical), triples, 
and quads with intervening Mgn absorption in one or all of the observed sight-lines. The mean redshift in our 
sample was (z a bs) = 1-099, covering a range of z = 0.043 — 2.066. 

Using a statistical approach, we tested whether the TC08 model can reproduce the observed variation in 
Mgn absorption strength between paired sight-lines. We varied the values of two model parameters: character- 
istic length £ c and covering fraction f c . As two sight-lines become arbitrarily close on the sky, we expect there 
to be some small separation below which the difference in Mgn absorption between sight-lines is negligible; 
we parameterize this effect as the characteristic length. The covering fraction measures the percentage of the 
gaseous halo projected on the sky that will produce Mgn absorption. 

The best fit parameters were £ c = 0.5 hy^ kpc and f c = 0.6 indicating that it is common to have large 
differences in magnitude of Mgll absorption between sight-lines that are separated by proper distances less 
than 0.5 /ifg 1 kpc. The covering fraction we find is inconsistent with that found by Chen & Tinker (2008) at 
80-86% if £ c > 2 h^} kpc. Our best-fit covering fraction is consistent with the 70% found by Chen et al. (2010) 
at z < 0.5, and the 50% found by Lundgren et al. (2011), though the latter measurement made with only one 
Mgn detection. 

We can rule out f c > 0.87 for all values of £ c at 99.7%. We can also rule out the region where £ c < 1.0 
kpc and f c < 0.3. 

This study used data from the literature only. There are a number of quadruply lensed quasars that do 
not have resolved spectra available for each individual line of sight. A short search in the literature finds 
SDSSJ1330+1810 (Oguri et al. 2008), SDSSJ1251+2935 (Kayo et al. 2007), WFI J2026-4536 and WFI J2033- 
4723 (Morgan et al. 2004). All four objects have either confirmed or suspected intervening Mgn absorption. 
It would be worth observing these objects spectroscopically with higher spatial resolution in order to separate 
out the individual lines of sight to see how the Mgn absorption varies between sight-lines. 

Ultimately, the largest uncertainty in this study is that we are ignorant of the properties and locations of 
the galaxies that are producing the intervening Mgn absorption seen in the spectra of quasar asterisms. A 
campaign to locate and study the galaxies responsible for absorption towards quasar asterisms would vastly 
improve models for Mgn absorption, such as the one studied here. 

Acknowledgements. We would like to thank Sara Ellison, Masamune Oguri, and Chris Churchill for providing 
data in support of this work. PBH and JAR acknowledge funding from NSERC, and JAR from the Faculty of 
Graduate Studies at York University. 

References 

[Barton & Cooke 2009] Barton, E. J. & Cooke, J. 2009, Astron. J., 138, 1817 



22 



[Bergeron & Boisse 1991] Bergeron, J. & Boisse, P. 1991, Astron. Astrophys., 243, 344 

[Bordoloi, Lilly, Knobel, Bolzonella, Kampczyk, Carollo, Iovino, Zucca, Contini, Kneib, Le Fevre, Mainieri, Renzini, Scodegg 
Bordoloi, R., Lilly, S. J., Knobel, O, Bolzonella, M., et al. 2011, arXi v:1106.06l"6l 

[Chen, Helsby, Gauthier, Shectman, Thompson, & Tinker 2010] Chen, H., Helsby, J. E., Gauthier, J., Sheet- 
man, S. A., Thompson, I. B., & Tinker, J. L. 2010, Astrophys. J., 714, 1521 

[Chen & Tinker 2008] Chen, H.-W. & Tinker, J. L. 2008, Astrophys. J., 687, 745 

[Churchill, Kacprzak, & Steidel 2005] Churchill, C. W., Kacprzak, G. G., & Steidel, C. C. 2005, in IAU Colloq. 
199: Probing Galaxies through Quasar Absorption Lines, ed. P. Williams, C.-G. Shu, & B. Menard, 24-41 

[Ellison, Ibata, Pettini, Lewis, Aracil, Petitjean, & Srianand 2004] Ellison, S. L., Ibata, R., Pettini, M., Lewis, 
G. P., Aracil, B., Petitjean, P., & Srianand, R. 2004, Astron. Astrophys., 414, 79 

[Kacprzak, Churchill, Barton, & Cooke 2011] Kacprzak, G. G., Churchill, C. W., Barton, E. J., & Cooke, J. 
2011, Astrophys. J., 733, 105 

[Kacprzak, Churchill, Steidel, & Murphy 2008] Kacprzak, G. C, Churchill, C. W., Steidel, C. C, k Murphy, 
M. T. 2008, Astron. J., 135, 922 

[Kayo, Inada, Oguri, Hall, Kochanek, Richards, Schneider, York, & Pan 2007] Kayo, I., Inada, N., Oguri, M., 
Hall, P. B., et al. 2007, Astron. J., 134, 1515 

[Lidman, Courbin, Kneib, Golse, Castander, & Soucail 2000] Lidman, C, Courbin, F., Kneib, J. -P., Golse, G., 
Castander, F., & Soucail, G. 2000, Astron. Astrophys., 364, L62 

[Lopez & Chen 2011] Lopez, G. & Chen, H.-W. 2011. larXiv:1110.3321l 

[Lopez, Hagen, & Reimers 2000] Lopez, S., Hagen, H.-J., & Reimers, D. 2000, Astron. Astrophys., 357, 37 
[Lovegrove & Simcoe 2011] Lovegrove, E. & Simcoe, R. A. 2011, Astrophys. J., 740, 30 

[Lundgren, Wake, Padmanabhan, Coil, & York 2011] Lundgren, B. F., Wake, D. A., Padmanabhan, N., Coil, 

A. , & York, D. G. 2011, Mon. Not. R. Astron. Soc, 417, 304 

[Mailer & Bullock 2004] Mailer, A. H. & Bullock, J. S. 2004, Mon. Not. R. Astron. Soc, 355, 694 

[Menard, Wild, Nestor, Quider, & Zibetti 2009] Menard, B., Wild, V., Nestor, D., Quider, A., & Zibetti, S. 
2009, Mon. Not. R. Astron. Soc, submitted. [arXiv:0912.3263l 

[Morgan, Caldwell, Schechter, Dressier, Egami, & Rix 2004] Morgan, N. D., Caldwell, J. A. R., Schechter, 
P. L., Dressier, A., Egami, E., & Rix, H. 2004, Astron. J., 127, 2617 

[Oguri, Inada, Blackburne, Shin, Kayo, Strauss, Schneider, & York 2008] Oguri, M., Inada, N., Blackburne, 
J. A., Shin, M., Kayo, I., Strauss, M. A., Schneider, D. P., & York, D. G. 2008, Mon. Not. R. Astron. Soc, 
391, 1973 

[Press, Teukolsky, Vetterling, & Flannery 2007] Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, 

B. P. 2007, Numerical Recipes 3rd Edition: The Art of Scientific Computing, 3rd edn. (New York, NY, USA: 
Cambridge University Press) 

[Rauch, Sargent, & Barlow 1999] Rauch, M., Sargent, W. L. W., & Barlow, T. A. 1999, Astrophys. J., 515, 
500 

[Spergel, Bean, Dore, Nolta, Bennett, Dunkley, Hinshaw, Jarosik, Komatsu, Page, Peiris, Verde, Halpern, Hill, Kogut, Limor 
Spergel, D. N., Bean, R., Dore, O., Nolta, M. R., et al. 2007, Astrophys. J. Suppl. Ser., 170, 377 

[Steidel, Dickinson, & Persson 1994] Steidel, C. C, Dickinson, M., & Persson, S. E. 1994, Astrophys. J. Lett., 
437, L75 

[Steidel, Erb, Shapley, Pettini, Reddy, Bogosavljevic, Rudie, & Rakic 2010] Steidel, C. C, Erb, D. K., Shap- 
ley, A. E., Pettini, M., Reddy, N., Bogosavljevic, M., Rudie, G. C, & Rakic, O. 2010, Astrophys. J., 717, 
289 



23 



[Tinker & Chen 2008] Tinker, J. L. & Chen, H. 2008, Astrophys. J., 679, 1218 

[Warren, Abazajian, Holz, & Teodoro 2006] Warren, M. S., Abazajian, K., Holz, D. E., & Teodoro, L. 2006, 
Astrophys. J., 646, 881 

[Zibetti, Menard, Nestor, Quider, Rao, & Turnshek 2007] Zibetti, S., Menard, B., Nestor, D. B., Quider, 
A. M., Rao, S. M., & Turnshek, D. A. 2007, Astrophys. J., 658, 161 



24 



